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Abstract 

We perform a detailed investigation of the scaling Potts field theory using the truncated 
conformal space approach. The aim is to check the scattering theory obtained from the S 
matrix bootstrap, both in the paramagnetic and ferromagnetic phases. It turns out that the 
bootstrap S matrices explain the finite size spectrum perfectly well. This finding is relevant 
in the view of a recent controversy surrounding the scaling limit of the lattice Potts model. 
In the conclusions we discuss possible scenarios to resolve the conflict between the perturbed 
CFT picture and lattice results, but the whole issue remains enigmatic. 

1 Introduction 

The three-state Potts model is one of the most important models in statistical mechanics. It is a 
generalization of the Ising model where the spin can take three values. At criticality, it is known 
to be integrable [IJ and the corresponding conformal field theory is also known P", '3]. 

Away from criticality, integrability on the lattice is lost. However in the vicinity of the critical 
point one can pass to a description by a continuum field theory, the scaling Potts model, which is 
uniquely specified by the operator content of the conformal field theory at the critical point [3] . 
The scaling Potts model is known to be integrable [4^ and its exact S matrix has been constructed 
using the bootstrap SI [SI [71 [5] , and also from an integrable lattice regularization or an integrable 
lattice regularization [3]. The bootstrap S matrix is diagonal i.e. there is only transmission without 
reflection. 

Recent investigations show that the scattering theory of the lattice model is not diagonal and 
it is suggested that it stays so even in the scaling limit [10]. A detailed numerical verification of 
this claim has also been carried out, and it turns out that the lattice spectrum does not agree with 
the predictions of the bootstrap S matrices |TT] (although there is partial agreement for a subset 
of the states). This is really surprising, since the identification of the scaling limit proceeds via 
well-known arguments that use the notion of universality of critical behavior. 

However, the bootstrap S matrices themselves have never been compared to the scaling Potts 
theory in detail. It is well-known that bootstrap S matrices are usually conjectures, and their 
correspondence to the microscopic (e.g. Lagrangian) description must be established. 

The thermodynamic Bethe Ansatz for the Potts model was constructed P^pl. but while it 
shows that the bootstrap scattering theory agrees in the ultraviolet with the conformal field theory 



detailed analysis of the scaling g-state Potts model scattering theory for general <? < 4, and the corresponding 
TBA system can be found in [T31IT1) . 
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predictions, it does not eventually test whether the S matrices really describe the scattering of 
multi-particle states. A simple way to see this is that the TBA equation is eventually the same 
for the Potts model and the tetra-critical Ising model, which is another perturbed conformal field 
theory with the same ultraviolet central charge, but the scattering in these two models show subtle 
differences, and they describe two different universality classes in terms of critical phenomena, with 
a different set of critical exponents. 

In view of this situation it is worthwhile to examine whether the bootstrap S matrices really 
describe the scattering theory of the scaling Potts model defined as a perturbed conformal field 
theory. This is the aim of this paper, accomplished using the truncated conformal space approach 
(TCSA) [12] to construct the finite volume spectrum of the scaling Potts field theory and comparing 
to the predictions made on the basis of the bootstrap S matrices. While the TCSA has been applied 
to the Potts model previously in [Tn|, that work mainly focused on the evolution of the particle 
spectrum in the non-integrable case when both the thermal and the magnetic perturbations are 
added. Here we stay with the purely thermal perturbation, but explore it much deeper in order 
to determine whether the conjectured bootstrap S matrices really describe the spectrum. 

The outline of the paper is as follows. In section 2 we collect some necessary facts about 
the scaling Potts model. Section 3 outlines the TCSA method, and in section 4 we discuss the 
predictions of the bootstrap S matrices for the finite volume spectrum. The detailed numerical 
comparison is made in section 5, while section 6 is reserved for the conclusions. 



2 Scaling Potts model as perturbed conformal field theory 

The three-state Potts quantum chain has three degrees of freedom per site. The Hilbert space is 
of the form 

i 

where i labels the sites of the chain, and the quantum space at site i has the basis |a) , a = 0, 1, 2. 
The dynamics of the Potts chain is defined by the Hamiltonian 

2 

i/chain - E E ^"^+1 - ^3 E (2.1) 

i a— i 

where 

P = I (i-'^oa' )!«)(«' I 

a.a' —0 

The critical point is at g = 1: for g > 1, the cyclic symmetry Z3 acting as 

\a) \a + 1 mod 3) 

is unbroken and the model is in a paramagnetic phase, while for 5 < 1 it is broken and the model 
is in a ferromagnetic phase. In the neutral sector of the model there is a duality under 

g^l/g (2.2) 

which is the analogue of the Kramers- Vannier duality in the Ising model. The critical point is 
self-dual and there the lattice model can be solved exactly [T]. 

The scaling (continuum) limit of the lattice Potts model is obtained by taking 

g -> 1 
such that (g - l)a"^/^ = fixed (2.3) 
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where g is the lattice Potts coupUng and a is the lattice spacing. Since the gap on the lattice 
scales as |10l [TT] 

A{g) « (g - 1)5/6 
this means a finite mass gap in the continuum limit 

a^^A{g) — > TO = finite 

On the other hand, if one sets the critical value 

5 = 1 

and takes the continuum limit a — > only afterward, then one arrives to a minimal conformal field 
theory with central charge 

4 



^=5 



The Kac table of conformal weights is 

{K,s} = 



2 1 1 4 7 ' 

\ 8 3 8 / 



s = 1, 



..,4 
..,5 



The D4 conformal field theory is invariant under the permutation group §3 generated by two 
elements Z and C with the relations 

= 1 = 1 czc = z-^ 

which have the signatures 

sign Z = +1 sign C = — 1 

The sectors of the Hilbert space are products of the irreducible representations of the left and right 
moving Virasoro algebras which can be specified by giving their left and right conformal weights 
as 

Sh.h = ® V,- 
The complete Hilbert space of the model is 



U = So.o ®S2 2 



551 j_ 

15 ' 15 



5 ' 5 

15 ' 15 



' ^3,3 



2 2 

3 ' 3 



352 7 

5 ' 5 



1 57 2 ® 5o,3 © 5; 



3,0 



(2.4) 



which corresponds to a modular invariant partition function on the torus [17J. Note that not 
all of the possible representations occur in the Hilbert space; there is another modular invariant 
partition function called A4 which includes all sectors of diagonal form Sh,h allowed by the Kac 
table exactly once: 



It corresponds to the scaling limit of a higher multicritical Ising class fixed point (with symmetry 
Z2). 

The sectors on the first line of (|2.4p are invariant under the action of the permutation group S3, 
the two pairs on the second line each form the two-dimensional irreducible representation, which 
is characterized by the following action of the generators: 



C\±) 
Z\±) 



IT) 



l±> 
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while the ones on the third hne form the one-dimensional signature representation where each 
element is represented by its signature. These sectors are in one-to-one correspondence with the 
families of conformal fields, and the primary field (the one with the lowest conformal weight) in 
the family corresponding to 5^ has left and right conformal weights h and h; they are denoted 
$^ with an optional upper ± index for fields forming a doublet of §3. In a family all other fields 
have conformal weights that differ from those of the primary by natural numbers. The conformal 
spin s = h — h gives the behavior under spacial translations; translational invariant fields must be 
spinless i.e. h = h. 

Relevant fields are exactly those for which h + h < 2. It is then obvious that the only §3- 
invariant spinless relevant field is 

$2 2 

which means that the Hamiltonian of the scaling limit of the off-critical Potts model is uniquely 
determined 

H = H^ + T J dx^2 2 (2.5) 
The coupling constant r is related to the lattice coupling in the scaling limit as 

T(x{g- gc)a'^'^ 

and r > in the paramagnetic, while r < in the ferromagnetic phase. The magnetization 
operator is the doublet field 



and there is also another relevant field doublet 

3 

which is the subleading magnetization operator. 

The scaling Potts model (j2.5p is integrable [4J, and its spectrum and the scattering matrix can 
be determined exactly [SJ HJ [SJ [71 IB] using the exact S matrix bootstrap developed in [TF] or an 
integrable lattice regularization [5] . In the paramagnetic phase, the vacuum is non-degenerate and 
the spectrum consists of a pair of particles A and A of mass m which form a doublet under §3 : 

c\A{e))^\A{e)) z\A{e))^e'r\A[e)) 

c\A{e)) = \A{e)) z\A{e)) = c-^|^(0)) (2.6) 

The generator C is identical to charge conjugation, the charge of A is Q = 1 while that of A is 
Q = —1 (where Q is defined mod 3). Due to two-dimensional Lorentz invariance, the energy and 
momentum of the particles is parametrized by the rapidity 6: 

E — m cosh 9 p = m sinh 9 

The two-particle scattering amplitudes are 



Saa{0i2) = SAAi012) 



sinh(^-l-f) 
sinh - f ) 



sinh + ^) 

Saa(.0,.) = ^^^(^-)--sinh(4-f) ^'-^^ 

where 9i2 — 9i—92 is the rapidity difference of the incoming particles. This S matrix was confirmed 
by the field theoretic thermodynamic Bethe Ansatz [12 . It can also be obtained by Bethe Ansatz 
from an integrable Z3 invariant lattice spin chain [5] , thought to be in the same universality class 
as the usual Potts chain defined by (|2.ip . 
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In the ferromagnetic phase the vacuum is three-fold degenerate 

|0)a a = -1,0,1 

where the action of §3 

Z\0)a=\0)a+lmod3 C\0)a = \0)^a 

and the excitations are kinks of mass m interpolating between adjacent vacua [3- The kink of 
rapidity 6, interpolating from a to 6 is denoted by 

Kabid) a-b = ±l mod 3 

The scattering processes of the kinks are of the form 

KabiOl) + KbM ^ KadiOl) + Kdc{02) 

with the scattering amplitudes equal to 

^ / \SAA{ei2) if a ^ ' 

The exact relation between the mass gap and the coupling r is also known |19| : 
|r| = Km^/^ 

_ r(^)[r(|)r(|)]«/^ /r(-i)r(i) 



4 X 2V5,s/^r (^) Vr(-|)r(|) 



0.1643033. 



3 TCSA for the scaling Potts model 

To test whether the S matrices in the previous subsection describe the scattering in the field 
theory correctly, we can turn to the truncated conformal space approach [15) . This consists of 
putting the system on a cylinder with compactified spatial dimension of circumference L and 
evaluating the excitation spectrum numerically. The matrix elements of the Hamiltonian between 
the r = eigenstates can be evaluated exactly using conformal field theory techniques. Continuing 
to Euclidean time Ie — —it and mapping the cylinder to the complex plane using 

z = eT(*^=+^^) , 2 = eT(*B-M) 

the explicit Hamiltonian is 
where is the matrix 

(i7o),, - (h, + - ^) (3.1) 
and B is given by the matrix elements of the perturbing field on the complex plane 

B,, = {i\^..{z^l)\j)5s,s, 

where Si and Sj are the conformal spins of the states \i) and \j) (the selection rule comes from 
translational invariance). The Hamiltonian can be made dimensionless in units of the mass gap 
m: 

1 27r / l^/^ \ 

4CSA = m^^^^^ = T (,^° ^ ''(2^^j 
/ — mL 
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where ± corresponds to the sign of r, i.e. to the para/ferromagnetic phases. To calculate the 
matrix elements B it is necessary to know the conformal three-point couplings which can be 
obtained from the literature [201 [5TJ [52] ; for the present calculations the construction in [53] was 
used. 

It turns out that the perturbing operator acts separately in the following four sectors: 

"Ho = So.O ® 52 2 © 57 2 © ^3.3 

5 ' 5 5 ' 5 

n± = s±^®st, 

15 ' 15 3 ' 3 

Hi = 52 7 ©57 2 ©5o,3©53,o (3.2) 

so the Hamiltonian can be diagonalized separately in each of them. It is also the case that the 
Hamiltonian is exactly identical in the sectors H+ and The reason for this is that the charge 
conjugation symmetry C acts on the sectors "H^ and "H- by swapping them. 

The spectrum is invariant under r — > — r in sectors T-Iq and T-Li- The latter fact is the conse- 
quence of a Z2 symmetry in these sectors under which the parities in T-Lq are 

even : 5o.o St 7 
5 ' 5 

odd : 52 2 53.3 

5 ' 5 

and so the perturbing operator $2 2 is odd. In "Hi this Z2 acts by swapping 

5o.3 -H- 53 52 7 o 5? 2 

' 5 ' 5 5 ' 5 

This symmetry leaves the fixed point Hamiltonian H^, and the conformal fusion rules in these 
sectors invariant^; away from the critical point, it is the continuum counterpart of the Kramers- 
Vannier duality (|2.2p under the scaling limit (|2.3p . 



4 Predictions for the finite volume spectrum 

4.1 Vacuum states in the ferromagnetic phase and vacuum tunnehng 

For T > one expects a single vacuum state, while for r < there should be three of them, 
corresponding to spontaneous breaking of the Z3 subgroup of §3. In a finite volume L, the 
degeneracies between these states must be lifted by tunneling. Using the formalism developed 
in |24| ■ tunneling in the three dimensional vacuum subspace can be described by the following 
effective Hamiltonian 

H=( 1 1 I (^— ^^e-™-^ + 0(e-2™-^) 
\ 1 1 / Vv27rmL 

where the matrix part is just the connectivity structure between the vacua introduced by the 
kinks, and the value of the matrix element is obtained by borrowing the dilute instanton gas result 
derived in [24J . The eigenvectors of this matrix are 



a=0,±l 
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with eigenvalues 



i?tu„nci(^9) = 2 cos 7? ^ e-"'^ + 0(e-2™^) 
v27rmL 



Two of the vacua remain degenerate 

^^tunncl(27r/3) — i?tunncl( — 27r/3) 



The conformal fusion rules do not allow the extension of this symmetry to the 'H± . 
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Figure 4.1: Vacuum adjacency diagrams for kinks in the ferromagnetic phase of the D4 and ^4 
models 

which is due to the presence of the charge conjugation Z2 subgroup of §3 which acts as 

C|±27r/3)L = |T27r/3)L 
C|0)l = |0)l 

therefore the degeneracy holds exactly (i.e. beyond the dilute instanton gas approximation). The 
splitting between these and the lowest lying state is given by 

A£;vac = £;tunnei(±27r/3) - E,^,,^,, (0) = 3 J^ e^"^ + 0(e~2™^) (4.1) 

v27rmL 

The vacuum states can also be described by thermodynamic Bethe Ansatz (TBA): the state |0)can 
be described by the TBA developed by Zaniolodchikov |25[ I12| . while the vacua | ± 27r/3) can be 
obtained by a twist [55]. The results from TBA also agree with the asymptotics (|4.ip . 

In connection with the TBA, note that the detailed structure of the A^ theory differs from 
that of Di also by including sectors absent from the latter. Although it has only one copy of the 
vacuum built upon 

15 ' 15 

it also contains a third vacuum state which in the ultraviolet limit corresponds to the field 

40 ' 40 

which is absent from the spectrum. Furthermore, while in the paramagnetic phase its scattering 
theory is identical to that of the Z?4, in the ferromagnetic phase its kink structure is markedly 
different [3] as shown in figure 14.11 

As a result of the above considerations, while TBA provides a very strong argument in favour 
of the correctness of the bootstrap S'-matrices, it cannot confirm their details since it cannot even 
distinguish between the ^4 and Z54 systems. Detailed confirmation is only possible by considering 
excited states, which is the motivation for the analysis that follows in the sequel. 



4.2 Multi-particle states in the paramagnetic phase 

The finite volume dependence of two-particle energy levels can be described using the standard 
Bethe- Yang equations, which are a specialization of Liischer's formalism for finite size corrections 
[771 to two-dimensional integrable field theories. The spectrum of multi-particle states 

\AaAei)---AaM) 

can be obtained by solving the Bethe- Yang equations which describe the periodicity of the wave 
function: 

g^n^Lsinhe, -Q^^^^^(^^ _^^.) ^ ^ fc=l,...n (4.2) 

The energy relative the vacuum state can be evaluated as 

n 

E{L)-Eo{L) = cosh 6lfc -l-0(e-''-^) 
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which is exact to all orders in 1/L i.e. up to corrections that vanish exponentially with the volume, 
as indicated. 

On the basis of Z quantum numbers, the states contained in sectors "Ho and "Hi are the Q = 
states 

AA , AAA , AAA , AAAA 

almost all of them occur in degenerate pairs under C; an exception is the state containing two 
stationary particles, A and A, as this is necessarily C-even. 

In finite volume, the C-even combinations are in "Hqj while the C-odd combinations in T-ix and 
the degeneracy, which is exact in the Bethe-Yang approximation, is lifted by corrections decaying 
exponentially with increasing volume. In addition, "Ho is expected to contain the vacuum, which 
is C-even. On the other hand, the states in H± are expected to be the Q = ±1 states 

: A,AA, AAA , AAAA 
U- : A, AA, AAA , AAAA 

and the two sectors are related to each other by C, with exactly identical spectra. 
4.3 Multi-particle states in the ferromagnetic phase 

In this case the kink structure necessitates the diagonalization of the multi-particle transfer matrix 
The two-particle states allowed by periodic boundary condition are of the form 

which is six dimensional corresponding to the allowed sequences 

aha = 010, 020, 121, 101, 202, 212 

and the states in the sector are overall neutral in terms of topological charge. The transfer matrices 
acting on this space, which describe the monodromy of the wave function corresponding to the 
kinks of rapidity d\ and Qi are 

Ti and 72 commute, and so they can be diagonalized simultaneously. The eigenvalues come in 
doubly degenerate pairs: 

g = +l,0,-l 

i.e. whenever Ti has eigenvalue t^^^ , 7i has eigenvalue tj'"* with the same eigenvector, and 
each of the three distinct eigenvalue pairs occurs twice corresponding to degeneracy under charge 
conjugation C. The quantization relations therefore read 

e'"'^'''"'^'^'e-'i^SAA{e2-ei) = 1 (4.3) 

This can be interpreted in the following way: there is a kink-particle correspondence 

Koi , K12 , K20 o A 
K02 , Kio , K21 ^ A 
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and there is a twist operator on the circle, which can be either the identity, Z or Z^^ . When 
calculating the monodromy of the wave-function, any kink taken around the circle gains the phase 
factor from the momentum and the scattering amplitudes with the other kinks which can be 
evaluated using the kink-particle correspondence. In addition there is a phase factor which is the 
eigenvalue of the corresponding particle under the twist operator, which accounts for the cases 
q = -1-1, 0, —1, respectively. 

This picture generalizes to higher particle states: the transfer matrix for the kink of rapidity 
Oi can be written as 

blb2...b^ _ Xb2Q fu^3 ^\ f,h 



ri(01,02, . . . ,en)ra..-. a„ = j (^1 " ^2)^ (^&3 «4 ) (^1 " 0,) 

X . . .b [ bn ai 

and due to the periodic boundary conditions the transfer matrices Tk for the other particles 
k = 2, . . . , n can be obtained by cyclically rotating the indices 1, . . . , n in this formula. All these 
transfer matrices commute and can be diagonalized simultaneously. An explicit computation shows 
that the three-particle states are described by 

e™^^"^'^e«"J'5AA(^i-^^2)5AA(0i-03) = 1 
e™^^"^'^e«^5AA(e2- 01)5^^(^2 -03) = 1 

e™^^"^'^e«^5AA(03-^^2)^AA(03-02) = 1 (4.4) 

and similar equations can be written for four-particle states. 

Due to periodic boundary conditions, all the states are of total charge zero. On the basis of 
their Z quantum numbers, the q = states are expected to be in Hq and Hi, coming in pairs under 
C which are degenerate in the Bethe-Yang approximation, so their degeneracy is only lifted by 
corrections exponentially decreasing with the volume. This is identical to the spectrum expected 
in the paramagnetic phase, in perfect agreement with the fact that the spectrum in Hq and Hi is 
invariant under the Kramers- Vannier duality r — >■ — r. 

The q = ±1 states are expected to be in the sectors 'H±, which again have exactly identical 
spectra due to the fact that C swaps them with other; these sectors also contain the vacua |±27r/3). 
As stated above, there are no states for which the total kink charge is nonzero; in particular, there 
are no one-kink states. Therefore these states are genuinely different from the paramagnetic 
case when sectors 'H± contained states with nonzero total charge. However, these sectors are 
characterized by the presence of the twist Z or Z~^, which in essence carries a twist quantum 
number q similar to the charge Q, distinguishing them from each other and the untwisted states. 



5 Analysis of the TCSA spectrum 

Here we show that the finite size spectrum derived from TCSA numerics is fully consistent with 
the theoretical picture presented in subsection |3] and so the detailed correspondence between the 
scattering matrices (|2.7l2.8p and the scaling field theory (j2.5p hold, both in the paramagnetic and 
the ferromagnetic regime. 

We diagonalized the Hamiltonian using the TCSA program developed for the calculations in 
[Jl]. It can be done separately in the four sectors. We cut off the Hilbert space by imposing an 
upper limit on the eigenvalues of the conformal Hamiltonian (j3.ip by considering only states with 

hi + hi< Ccut 

and restricted to the zero-momentum sector in which 

hi — hi = 

The dimensions of the Hilbert spaces for the values of the cutoff used in the calculation are 
summarized in the following table: 
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Figure 5.1: The gap in the ferromagnetic and paramagnetic phases. The dots are the TCSA data, 
while the continuous hne shows the prediction (|4.ip for the vacuum spht and (|5.ip for the volume 
dependent particle mass, respectively. 





16 


18 


20 


22 


24 


26 


Ho 


395 


744 


1438 


2591 


4708 


8174 


■Hi 


384 


732 


1422 


2574 


4682 


8144 


H± 


600 


1225 


2262 


4287 


7652 


N/A 



One technical issue is that the TCSA converges slowly; however, it is ultraviolet finite due to 
the fact that the conformal weight of the perturbing operator is less than 1/2. Convergence can 
be improved using the TCSA renormalization group [32i [551 which predicts that the cutoff 
dependence of energy levels computed from TCSA is 

Ei{L, Bcut) = Ei{L, oo) + bi{L) h subleading corrections 

Using this result, the computed spectra can be extrapolated to infinite cutoff, which results in a 
substantial improvement of accuracy. 



5.1 The gap 



First we analyze the gap in the spectrum, which is the difference between the ground state energy 
in sectors 'H± and Hq. As shown in figure 15.11 the gap vanishes exponentially with L in the 
ferromagnetic phase in accordance with (j4.ip . which corresponds to a triply degenerate vacuum 
related to spontaneously broken Z3 symmetry. In the paramagnetic phase the gap approaches 
unity, i.e. in proper units m as L grows, corresponding to the one-particle states A and A, which 
are expected to be the lowest lying states in sectors 'H±. For a more precise comparison, one can 
also use Liischer's formalism for finite volume mass corrections [35] (see also [55]) which gives the 
following theoretical prediction for the finite volume mass gap in the paramagnetic phase: 



AE{L) 



TO - /ic/^e '^^ 




Saa 



AA 



X m cosh Oe 



mlj cos 



he 



O 



(5.1) 



Res Saa{0) 

=27ri/3 



V3 



x/3 



The numerical accuracy of the agreement with the expected gap in the range 2.5 < mL < 10 is of 
order 10"'^ or better. This gives a good test of the S matrices, but again, it is not really sensitive 
to details such as whether one deals with the D4 or the A4 model. It also does not give a test of 
the individual two-particle iS-matrices separately. 
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AE{L)/m 




mL 



Figure 5.2: The spectrum in the Hq sector 



5.2 Excited states in "Hq ^ind Hi 

To test the S matrices, we first examine the sector Hq. Figure shows the agreement between 
the two-particle hne predicted by the Bethe-Yang equations from the bootstrap S matrices and 
the TCSA data. The continuous hues are two-particle states, but to explain all the states shown 
it is also necessary to consider three-particle AAA levels (dashed lines) and four-particle AAAA 
levels (dotted lines), which could be done by a simple extension of the formalism formulated 
in subsection S) The lowest lying two-particle level contains stationary A and A particles; the 
Bethe-Yang prediction for this level is constant: 

AE{L) = 2m 

but it has exponential corrections similar to those for the mass gap in (I5.ip . which clearly show 
up on the plot. It only occurs in the C-even sector Hq. 

The agreement is spectacular, but we must keep in mind that this sector is independent of the 
phase, and even of the choice of the universality class A4 or D4. 

The same agreement is observed in sector Hi shown in figure [531 but this is more interesting, 
given that this sector occurs only in the D4 theory. Comparing the numerical spectra it can also 
be seen that states with the same particle contents are almost degenerate between Hq and "Hi, 
split only by corrections decaying exponentially with L. 

5.3 n± for r > 

In the paramagnetic phase, the two-particle and three-particle levels can be interpreted in terms 
of AAA and AAA states in the 7i+ and their charge conjugates in the H- sector, as expected 
from the description given in subsection 01 This is demonstrated in figure [531 which shows one of 
the sectors, as the spectrum in the other one is identical. 

5.4 n± for r < 

In the ferromagnetic phase, just as expected, this sector turns out to contain the two-kink states 
described by the g = ±1 case of the Bethe-Yang equations (|4.3p . This is shown in figure [5751 
where the two-particle lines are shown by continuous lines. The other states in the plot can 
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AE(L)/m 




niL 



Figure 5.3: The spectrum in the T-Li sector 



AE(L)/m 



20 r 



15 - 



10 - 



5 - 




mL 



Figure 5.4: The spectrum in the T-L± sectors for T > Tc- The lowest lying (one-particle) level has 
been omitted from the plot, since it was already analyzed in detail in 15. II 
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AEitym 



20 r 




15 



10 



10 



mL 



Figure 5.5: The spectrum in the H± sectors for T < Tc- The lowest lying (vacuum) level has been 
omitted from the plot, since it was already analyzed in detail in 15. II 

be explained using the three-kink equation (j4.4l) and the four-kink generalization of the transfer 
matrix formalism, and are shown by dashed and dotted lines, respectively. 

5.5 Explicit determination of the phase-shift functions 

Although the agreement between the Bethe-Yang predictions and the TCSA spectrum seems quite 
convincing, it is not obvious whether this is true for the phase-shifts, since the dependence on the 
phase-shift is a subleading term in the finite size corrections. However, we can extract the phase- 
shifts directly from two-particle states. For AA we can consider the states in the zero-momentum 
subspaces of "Ho and "Hi (identical in both regimes) which are described by 



that follows from taking the logarithm of the appropriate Bethe-Yang equation. The phase-shift 
function can be extracted using 



and plotted as a parametric curve in the variable / — mL. The two sectors give exponentially close 
results, which are visually indistinguishable, so only one of them is shown in figure [5TB1 (a). It is 
also possible to extract the same phase shift function from the paramagnetic 'H± sectors, which 
are described by the same equations, but with 



The numerical results from this extraction are shown in figure 15.61 (b) . 

For the case of AA, only the paramagnetic phase T-L± sectors have such two-particle states 
because these are charged. The phase-shifts determined from the first two such levels are shown 
in figure 15.71 For large values of 9 (corresponding to small mL) the observed deviations can be 



mL sinh6' -I- 5aa{'^(^) 



27r/ / e N 
-tlogS^AiO) 



9{L) 
SAAmL)) 



cosh~^ {AE{L)/2m) 
27rl ~ mLsmhe{L) 
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S{0) 6(8) 




(a) (b) 

Figure 5.6: The phase-shift function Sj^^{d) = —i log 5*^^(6*) determined from the first two zero- 
momentum two-particle levels in sectors (a) Hq and (b) in the ferromagnetic regime. The 
continuous lines are the theoretical prediction, while the dots show the values obtained from 
TCSA. 



S(6) 




Figure 5.7: The phase-shift function 6aa{S) = —i\ogSAA{0) determined from the first two zero- 
momentum two-particle levels in the sector in the paramagnetic regime. The continuous lines 
are the theoretical prediction, while the dots show the values obtained from TCSA. 



attributed to finite-size effects decaying exponentially with the volume (which are neglected in the 
Bethe-Yang description) , while for small values of 6 (large mL) they reflect the truncation error 
inherent in the TCSA method. 

Note that these phase-shifts enter in the determination of many other levels, such as three and 
four-particle states, for which the agreement between the theoretical prediction and TCSA data 
is as good for the two-particle levels. This further strengthens the agreement, and also confirms 
that the AA phase-shift correctly describes the scattering of like-charge kinks in the ferromagnetic 
phase, although a direct determination from two-particle states (due to their absence from the 
T < Tc spectrum) is not possible. 

6 Conclusions and open questions 

The main conclusion is very simple: the S matrices (|2.7I) and (12. 8p do describe the scattering 
theory of the r > and r < phase of the scaling Potts model, defined as the perturbed 
conformal field theory (|2.5p . Everyone familiar with the field of integrability would say that this 
was to be expected. However, this leaves us with an enigma: why is it that these S matrices are 
incompatible with the lattice results as shown in [TUl 111) ? 

The simplest scenario is that the lattice calculations are simply too far from the critical point, 
and so not in the scaling regime yet. The results in ^llj however seem to indicate that that the 
scaling regime of the lattice model was reached at least for the low-energy degrees of freedom; the 
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lattice calculation also reproduces the correct S'-matrix for the C-odd neutral states (sector 
Hi). The lattice data however show no trace of the exponential degeneracy between the C-even 
multi-particle states in sector Hq and their C-odd counterparts in sector "Hi, and yield a scattering 
length for the C-even states that is inconsistent with Sj^^. In this scenario this could be due 
to the contribution of some irrelevant operator with an unusually large coefRcient, and when 
approaching the critical point on the lattice the behaviour would be described by the scaling Potts 
model for a range of energy values which excludes very small rapidities (and also large ones, due 
to the presence of the cut-off), but extends longer and longer as 5 — > 1 such that the scattering 
amplitude determined from the lattice tends to everywhere (except at ^ = 0, where it always 
remains —1, whereas >S'^^(0) = +1). 

The above scenario is, however, not really supported by the lattice data in pT] (although it 
is not ruled out conclusively, either). Therefore there is also another, more exotic scenario: it is 
possible that the low energy behaviour of the off-critical lattice Potts model is not the integrable 
scaling Potts model, but some other (probably non-integrable) field theory X. We know that the 
critical Potts model in the scaling limit is described by the c = 4/5 CFT, which has a single 
Sa-invariant relevant operator, so X cannot have the c = 4/5 D4 CFT as its ultraviolet limit. If 
this scenario makes any sense, the open question is: what is theory X? 
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